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Abstract. Dust ion acoustic solitary structures have been investigated in an 
unmagnetized nonthermal plasma consisting of negatively charged dust grains, 
adiabatic positive ions and nonthermal electrons. For isothermal electrons, the 
present plasma system does not support any double layer solution, whereas for 
nonthermal electrons, negative potential double layer starts to occur whenever the 
nonthermal parameter exceeds a critical value. However this double layer solution 
is unable to restrict the occurrence of all negative potential solitary waves of the 
present system. As a result, two different types of negative potential solitary waves 
have been observed, in which occurrence of first type of solitary wave is restricted 
by Mc < M < Mjj whereas the second type solitary wave exists for all M > M/j, 
where Mc is the lower bound of Mach number M, i.e., solitary structures start 
to exist for M > Mc and Mo(> Mc) is the Mach number corresponding to a 
negative potential double layer. A finite jump between the amplitudes of negative 
potential of solitary waves at M = Mo — ei and at M = Mo + 62 has been 
observed, where < ei < Mo — Mc and £2 > 0. As double layer solution plays 
an important role for the present system, an analytical theory for the existence 
of double layer has been presented. A numerical scheme has also been provided 
to find the value of Mach number at which double layer solution exists and also 
the amplitude of that double layer. The solitary structures of both polarities can 
coexist whenever fi exceeds a critical value, where fi is the ratio of the unperturbed 
number density of electrons to that of ions. Although the occurrence of coexistence 
of solitary structures of both polarities is restricted by Mc < M < Mmax , only 
negative potential solitary wave still exists for all M > Mmax, where M^ax is 
the upper bound of M for the existence of positive potential solitary waves only. 
Qualitatively different solution spaces, i.e., the compositional parameter spaces 
showing the nature of existing solitary structures of the energy integral have been 
found. These solution spaces are capable of producing new results and physical 
ideas for the formation of solitary structures whenever one can move the solution 
spaces through the family of curves parallel to the curve M = Mc- 
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1. Introduction 

Acoustic wave modes in dusty plasma have received a great deal of attention since 
the last decade [TJ[2l[3l|4l[5l|6l[7l[8]. Depending on different time scales, there can 
exists two or more acoustic waves in a typical dusty plasma. Dust Acoustic (DA) and 
Dust Ion- Acoustic (DIA) waves are two such acoustic waves in a plasma containing 
electrons, ions, and charged dust grains. 

Shukla and Silin [5] were the first to show that due to the quasi neutrality 
condition n^Q + tliiqZ^ = riio and the strong inequality n^o <C n^o (rieo, n^o, and 
Udo are, respectively, the number density of electrons, ions, and dust particles, where 
Zd is the number of electrons residing on the dust grain surface), a dusty plasma (with 
negatively charged static dust grains) supports low-frequency DIA waves with phase 
velocity much smaller (larger) than electron (ion) thermal velocity. In case of long 
wavelength limit the dispersion relation of DIA wave is similar to that of Ion- Acoustic 
(lA) wave for a plasma with Ueo = nm and Ti -^Te, where Ti(Te) is the average ion 
(electron) temperature. Due to the usual dusty plasma approximations (rieo ^ ?^^o 
and Ti c:i T^), a dusty plasma cannot support the usual I A waves, but the DIA waves 
of Shukla and Silin [2j can. Thus DIA waves are basically lA waves, modified by the 
presence of heavy dust particulates. The theoretical prediction of Shukla and Silin [2] 
was supported by a number of laboratory experiments [^ 1101 [TT) . The linear properties 
of DIA waves in dusty plasma are now well understood [HI [T31 [Tl] . 

Dust Ion-Acoustic solitary waves (DIASWs) have been investigated by several 
authors. Bharuthram and Shukla [15] studied the DIASWs in an unmagnetized 
dusty plasma consisting of isothermal electrons, cold ions, in both static and mobile 
dust particles. Employing reductive perturbation method, Mamun and Shukla [l6] 
investigated the cylindrical and spherical DIASWs in an unmagnetized dusty plasma 
consisting of inertial ions, isothermal electrons, and stationary dust particles. They 
[17] have also investigated the condition for existence of positive and negative potential 
DIASWs. Verheest et al. [18] have shown that in the dust-modified ion acoustic 
regime, negative structures can also be generated, beside positive potential soliton 
if the polytropic index 7e 7^ 1 for electrons. The effect of ion-fluid temperature on 
DIASWs structures have been investigated by Sayed and Mamun in a dusty plasma 
containing adiabatic ion-fiuid, Boltzmann electrons, and static dust particles. 

In most of the earlier works, Maxwellian velocity distribution function for lighter 
species of particles has been used to study DIASWs and DIA double layers (DIADLs) . 
However the dusty plasma with nonthermally /suprathermally distributed electrons 
observed in a number of heliospheric environments [T31 dU [201 1211 1221 123] • Therefore, 
it is of considerable importance to study nonlinear wave structures in a dusty plasma in 
which lighter species (electrons) is nonthermally /suprathermally distributed. Bcrbri 
and Tribeche [24] have investigated weakly nonlinear DIA shock waves in a dusty 
plasma with nonthermal electrons. Recently Baluku et al. |25j have investigated 
DIASWs in an unmagnetized dusty plasma consisting of cold dust particles and kappa 
distributed electrons using both small and arbitrary amplitude techniques. 

In the present investigation we have considered the problem of existence of 
DIASWs and DIADLs in a plasma consisting of negatively charged dust grains, 
adiabatic positive ions and nonthermal electrons. Three basic parameters of the 
present dusty plasma system are /i, a and which are respectively the ratio of 
unperturbed number density of nonthermal electrons to that of ions, the ratio of 
average temperature of ions to that of nonthermal electrons, a parameter associated 
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with the nonthermal distribution of electrons. Nonthermal distribution of electrons 
becomes isothermal one if (3i = 0. 

The main aim of this paper is to investigated DIASWs and DIADLs thoroughly, 
giving special emphasis on the followings: 

(a) To study the nonlinear properties of DIA waves in a dusty plasma with nonthermal 
electrons, (b) To find the exact bounds (lower and upper) of the Mach number AI 
for the existence of solitary wave solutions, (c) As double layer solution plays an 
important role to restrict the occurrence of at least one sequence of solitary waves 
of same polarity, we set up an analytical theory to find the double layer solution 
of the energy integral, which help us to find the Mach number at which double 
layer occurs and also, to find the amplitude of that double layer solution, (d) On 
the basis of the analytical theory for the existence of solitary waves and double 
layers, the present plasma system has been analyzed numerically. Actually, analyzing 
the Sagdeev potential, we have found qualitatively different solution spaces or the 
compositional parameter spaces showing the nature of existing solitary structures 
of the energy integral. From these solution spaces, the main observations are the 
followings. (dl) For isothermal electrons, the present plasma system does not support 
any double layer solution in both cold and adiabatic cases. For nonthermal electrons, 
the present plasma system does not support any Positive Potential Double Layer 
(PPDL) solution, whereas Negative Potential Double Layers (NPDLs) start to occur 
whenever the nonthermal parameter exceeds a critical value. However this NPDL 
solution is unable to restrict the occurrence of all Negative Potential Solitary Waves 
(NPSWs) of the present system, i.e., NPDL solution is not the ultimate solution of 
the energy integral in the negative potential side. Actually, we have observed two 
different types of NPSWs, in which amplitude of first type of NPSW is restricted 
by the amplitude of NPDL whereas the amplitude of NPDL is unable to restrict the 
amplitude of the second type NPSWs. As a result, we have observed a finite jump in 
amplitudes between two different types of NPSWs separated by a NPDL. This fact 
has also been observed recently by Verheest[26j and Baluku et al.^7\ for ion-acoustic 
solitary wave with different plasma constituents. (d2) For any physically admissible 
values of the parameters of the system, specifically, for any value of and any value 
of (3i, NPSW exists for all M > Mc except M = Md, where is the lower bound of 
Mach number M, i.e., solitary structures start to exist for M > Mc and Md{> Mc) 
is the Mach number corresponding to a NPDL solution. However, if the parameter 
ji exceeds a critical value ^p, PPSWs exist for all < ^ < ^t(< 1) whenever the 
Mach number lies within the interval Mc < M < M^ax, where /^r(< 1) is a physically 
admissible upper bound of /z and M^ax is the upper bound of M for the existence 
of PPSWs only, i.e., there does not exist any PPSW if M > M^ax- Therefore, the 
coexistence of both PPSWs and NPSWs is possible for all /ip < /i < /ir whenever 
Mc < M < Mmax, but NPSWs still exist for M > M^ax- (d3) For nonthermal 
electron species, we have investigated the entire solution space of the energy integral 
with respect to the nonthermal parameter /3i and we have found four qualitatively 
different solution spaces depending on the cut off values of ^. Actually, here we are 
able to define three cut off values /j,p, /Xg and of n such that Q < iip < < ^.r < 1 
for any given value of a, and consequently, we can partition the entire range of /i in 
the following four disjoint subintervals: < /i < /Xp, /ip < /i < /ig, /ig < < Mr ^-nd 
Mr ^ M ^ Mt(< !)• For these four disjoint subintervals of /x, we have four different 
solution spaces of the energy integral with respect to nonthermal parameter j3i . These 
solution spaces can define all types of solitary structures of the present system. (d4) 
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Finally, considering any solution space, one can get new results and physical ideas 
for the formation of solitary structures if he moves the solution space through the 
family of curves parallel to the curve M = M^- If we move the solution space through 
the family of curves parallel to the curve M = Mc, it is simple to understand the 
mathematics as well as physics for the formation of double layer solution and it is also 
simple to understand the relation between solitons and double layer solution. 

The present paper is organized as follows: Basic equations are given in section [2] 
Derivation of energy integral along with Sagdeev potential is given in section [3l 
Physical interpretation for the existence of solitary structures of the energy integral is 
given in section |4l The lower and upper bounds of the Mach number for the existence 
of solitary structures are given in section [5l In section 15.11 the analytical method to 
find the upper bound of the Mach number for the existence of PPSWs is given. In 
section 15.21 we find that the existence of NPDL solution may restrict occurrence 
of NPSWs having amplitude less than the amplitude of NPDL. In section 15.31 
an analytical theory to find the double layer solution of the energy integral has 
been provided. In section 15.41 an algorithm has been provided to find the value 
of Mach number at which double layer solution exists and also the amplitude of 
that double layer. Following a logical sequence of numerical scheme based on the 
theoretical discussions as given in section |4] and section [SJ the solution spaces have 
been constructed in section [5] Finally, we have concluded our findings in section [71 

2. Basic equations 

The governing equations describing the nonlinear behavior of DIA waves, propagating 
along X-axis, in collisionless, unniagnetized dusty plasma consisting of negatively 
charged immobile dust grains are the following: 

rii^t + {riiUi)^ = 0, (1) 

Uij. + UiUi^x = -(t>x - an~^pi^x, (2) 

Pi,t + UiPi^x + lPiUi,x = 0, (3) 

<Pxx = {neo/nio)ne - rii + (^d'^do)/»^zO, (4) 
where the parameter a — Ti/T,,. 

Here we have used the notation V'i, q or {tl'ijq for dipi/dq and n^, rig, Ui, pi, 
cj), X and t are, respectively, the ion number density, electron number density, ion 
velocity, ion pressure, electrostatic potential, spatial variable and time, and they have 
been normalized by n^o (unperturbed ion number density), rieo (unperturbed electron 
number density), Ci{— \J {K bT(,) / rrii) (ion-acoustic speed), UioKsTi, KsTe/e, 
^Dem{= a/ (ii'Bre)/(47rnioe2))(Debye length), and uj~^{= rrn / {Aiiniae^)) (ion 
plasma period). Here 3) is the adiabatic index, Kb is the Boltzmann constant, 
Ti and Tg are, respectively, the average temperatures of ions and electrons, is the 
mass of an ion, n^o is the dust number density, Zd is the number of negative unit 
charges residing on dust grain surface, and e is the charge of an electron. 

The above equations are supplemented by nonthermally distributed electrons as 
prescribed by Cairns et al [28l for the electron species. Nonthermal distribution of any 
lighter species of particles (as prescribed by Cairns et al [28] for the electron species) 
can be regarded as a modified Boltzmannian distribution, which has the property 
that the number of particles in phase space in the neighborhood of the point v = Q 
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is much smaller than the number of particles in phase space in the neighborhood of 
the point t; = for the case of Boltzmann distribution, where v is the velocity of the 
particle in phase space. Under the above mentioned normalization of the dependent 
and independent variables, the normalized number density of nonthermal electrons 
can be written as 

n, = (1 ^ M + M'')et (5) 

where 

with ai > 0. Here /3i{ai) is the parameter associated with nonthermal distribution 
of electrons and this parameter determines the proportion of fast energetic electrons. 
From ^ and the inequality ai > 0, it can be easily checked that the nonthermal 
parameter /3i is restricted by the following inequality: < /3i < 4/3. However 
we cannot take the whole region of /3i (0 < /3i < 4/3). Plotting the nonthermal 
velocity distribution of electrons against its velocity (v) in phase space, it can be 
easily shown that the number of electrons in phase space in the neighborhood of the 
point V = decreases with increasing /3i and the number of electrons in phase space 
in the neighborhood of the point w = is almost zero when /3i — )■ 4/3. Therefore, 
for increasing /3i distribution function develops wings, which become stronger as /3i 
increases, and at the same time the center density in phase space drops, the latter 
as a result of the normalization of the area under the integral. Consequently, we 
should not take values of /?i > 4/7 since that stage might stretch the credibility of the 
Cairns model too far [53]. So, here we consider the effective range of /3i as follows: 
< /3i < Pit, where /3it = 4/7 « 0.571429. 
The charge neutrality condition, 

Zdrido + n-eo = Uio, (7) 

can be written as 

^ = 1-M, (8) 
riiQ 

where ^ = rieo / and consequently, the Poisson equation (|4]) assumes the following 
form: 

(l)xx = IJ-ne-rii + il- fi). (9) 

We note from ([8]) that 1 — /i must be greater than zero, i.e., < ^ < 1. When 
/Li — >■ 1, the effect of negatively charged dust grains on DIA wave is negligible and so, 
we restrict fi by the inequality < /i < fix, where fj,T is strictly less than 1. 

3. Energy integral and the Sagdeev potential 

To study the arbitrary amplitude time independent DIASWs and DIADLs we make 
all the dependent variables depend only on a single variable — x — Mt, where the 
Mach number M is normalized by c^. Thus in the steady state, (HI - © and ^ can 
be written as 

- Mn,, 5 + {fiiu,)^ = 0, (10) 

- Mu,^ 5 + u,Ui^ ^ = - an~^pi^ (11) 
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- Mpi^ ^ + UiPi^ ^ + ipiUi^ 5 = 0, (12) 

(f)^^ = fine - + (1 - fi). (13) 
Using the boundary conditions, 

Hi I, Pi -> 0, 

(/) 0, 0^ ^ as 1^1 oo, (14) 

and solving (|10p - (|12p . we get a quadratic equation for nf and following the same 
argument as given in Das et al. |30j to find the expression of dust density with exact 
bounds, we get the following expression of n^. 

y/2M 



(15) 



where 

Wa/ ^ • (Id) 

From (jl5l) . we see that this equation gives both theoretically and numerically correct 
expression of even when a = if < '5 a/ • 

Now integrating (fT3)) with respect to and using the boundary conditions ([14]), we 
get the following energy integral with V{(j)) as Sagdeev potential or pseudo-potential. 

where 

!/(,/)) = V{M, cj))^V,^ fiVe - (1 - fi)cb, (18) 
- (1 + 3/3i - 3/3i0 + /3i02)e"^ - (1 + 3/3i), (19) 

+ a- ni{M^ + 3a - 2(j) - 2anf). (20) 

Here, t^(M, (/>) is same as V^(0), i.e., the Mach number M is omitted from the notation 
y(M, (/)) when no particular emphasis is put upon it. 



4. Physical interpretation of the energy integral 

The energy integral (jl7p can be regarded as the one-dimensional motion of a particle 
of unit mass whose position is at time ^ with velocity dcf)/ dS^ in a potential well V{(j)). 
The first term in the energy integral (IT71) can be regarded as the kinetic energy of a 
particle of unit mass at position (j) and time ^ whereas V{(j)) is the potential energy 
at that instant. Since kinetic energy is always a non-negative quantity, V{(j)) < for 
the entire motion, i.e., zero is the maximum value for V{(j)). Again from (jl7p . we find 
d'^(j>/d^'^ ~\- V'{(j)) — 0, i.e. the force acting on the particle of unit mass at the position 
(f) is — y (0), where " ' " indicates a derivative with respect to (j). Now, it can be easily 
checked that V^(0) = V'{Q) — 0, and consequently, the particle is in equilibrium at 
= because the velocity as well as the force acting on the particle at = are 
simultaneously zero. Now if = can be made an unstable position of equilibrium, the 
energy integral can be interpreted as the motion of an oscillatory particle \iV{(l)m) = 
for some 7^ 0, i.e., if the particle is slightly displaced from its unstable position 
of equilibrium then it moves away from its unstable position of equilibrium and it 
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continues its motion until its velocity is equal to zero, i.e., until cf) takes the value 0^. 
Now the force acting on the particle of unit mass at position = 0,„ is —V'{(f>m). For 
4>m < 0, the force acting on the particle at the point = 0^ is directed towards the 
point = if —V'{(f>m) > 0, i.e., if V'{(j)m) < 0. On the other hand, for 0™ > 0, 
the force acting on the particle at the point (j) — 0™ is directed towards the point 
= if -V'((l)m) < 0, i.e., if > 0. Therefore, if V'{(j)m) > (for the positive 

potential side ) or if V'{(j)m) < (for the negative potential side ) then the particle 
reflects back again to = 0. Again, if F(0„i) = V'{(l)m) — then the velocity dcj)/d^ 
as well as the force (Pcfi/d^,'^ both are equal to zero at (j) — (j)„i. Consequently, if 
the particle is slightly displaced from its unstable position of equilibrium (0 — 0) 
it moves away from = and it continues its motion until the velocity is equal 
to zero, i.e., until takes the value = 0^. However it cannot be reflected back 
again at = as the velocity and the force acting on the particle at = 0,„ vanish 
simultaneously. Actually, if T^'(0m) > (for 0™ > 0) or if T^'(0m) < (for 0^ < 0) 
the particle takes an infinite long time to move away from the unstable position 
of equilibrium. After that it continues its motion until takes the value 0m and 
again it takes an infinite long time to come back its unstable position of equilibrium. 
Therefore, for the existence of a positive (negative) potential solitary wave solution 
of the energy integral (|17|) . we must have the following: (a) = is the position 
of unstable equilibrium of the particle, (b) y(0m) = 0, F'(0,„) > {V'{(j)m) < 0) 
for some 0™ > (0m < 0), which is nothing but the condition for oscillation of the 
particle within the interval min{0, 0^,} < < max{0, 0m} and (c) T^(0) < for all 
< < 0m (0m < < 0), which is the condition to define the energy integral (fTTl) 
within the interval min{0, 0m} < < max{0, 0m}. For the existence of a positive 
(negative) potential DL solution of the energy integral pT|) . the conditions (a) and 
(c) remain unchanged but here (b) has been modified in such a way that the particle 
cannot be refiected again at = 0, i.e., the condition (b) assumes the following form: 

Vi(l)m) = V'icj,m) = 0, y"(0m) < for SOmC 0m > (0m < 0). 



5. Lower and Upper bounds of Mach number for the existence of solitary 
waves and double layers 

For the existence of solitary structures, we must have V{0) = V'{0) — and V"{0) < 0. 
Now it can be easily verified that the first two conditions, i.e., V^(0) = and V'{0) = 
are trivially satisfied whereas, the condition V"{0) < gives M > M^, where is 
given by the following equation. 

Mc=3a+— 1— . (21) 

Now, for Mc to be real and positive, we must have > and < /3i < 1. As the 
effective range of ;3i is < /3i < /3it, where (3it = 4/7 « 0.571429, is well-defined 
as a real positive quantity for all < /i < /zt and < /3i < /3it- 



5.1. Upper bounds of the Mach number for the existence of PPSWs 

Consider the existence of a PPSW for some value of AI > Ale. Therefore, there exists 
a 0m > such that 



1/(0) < for all < < 0m, 

Vi^m) = , ■l/'(0m) > 0. 



(22) 
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Now as is real for cj) < \l/jv/ we must have (j)m < ^'m, otherwise V{(j)jn) is not a 

real quantity. Therefore, 

V{(t)) < for all 0<4>< (0„ <)*M, (23) 

defines a large amplitude PPSW, which is in conformity with (|22]) . if F(^'m) = and 
V'i^M) > 0. 

Again let M„iax be the maximum value of M up to which solitary wave solution 
can exist. As '^'m increases with M then '^'m < ^m™„j; ■ Therefore, 

< for aU < 0< < *m <)*m,„„., 

defines the largest amplitude PPSW if F(*a/„„J = and y'(*A/„„J > 0. Therefore, 
for the existence of the PPSWs, the Mach number M is restricted by the following 
inequality: < AI < M^ax^ where M^ax is the largest positive root of the equation 
V{'^m) = subject to the condition T^(*m) > for aU M < M^ax- Now if at 
M = Md (Mc < Md < Mjnax), one can get a PPDL solution of the energy integral 
P7I) then for the existence of the PPSWs, the Mach number AI is restricted by the 
inequality: Afc < Ad < Add and also Aid < AI < Almax- On the other hand if 
Aid — M,nax, then for the existence of PPSWs, the Mach number M is restricted by 
the inequality: Ale < AI < Md- 

5.2. Upper hounds of the Mach number for the existence of NPSWs 

We have seen earlier that V{4)) is real if < m, where m is strictly positive. For 
NPSWs or NPDLs, we have 

V{(f)) < for aU < < 0, (24) 

along with the other conditions stated in section U for the existence of NPSWs or 
NPDLs. As is strictly positive and for NPSWs or NPDLs, (j) < 0, the condition 
< \1/m is automatically satisfied and consequently, for these two cases V{4>) is 
well defined for all < without imposing extra condition. Since there is no such 
restriction on 0, we cannot use the same definition as in the case of PPSWs to find 
the upper bound of Mach numbers for the existence of NPSWs. For the case of 
NPSWs, to find an upper limit or upper bound of AI, up to which NPSW can exist, 
we shall first of all find a value AIjj of AI for which energy integral p7|) gives a NPDL 
solution at AI — AIo with amplitude (p — (po- Now if at AI = Air,, (t> = <t>D the 
only root (double root) of the equation V{<j)) = y(Af, 0) = 0, i.e., V{AIdt4'd) — 
and V'{Mm <t>D) = 0, then the NPDL solution is the ultimate solution of the energy 
integral ([TT]) and in this case, no NPSW solution can be obtained if M > AIo, i-e., 
for the occurrence of NPSWs, the Mach number M is restricted by the inequality 
Mc < M < AId- 

On the other hand if there exists an inaccessible simple root 0c i of such that 
(pDi < (t)D, i.e., V{AID,(t)D^) = along with T/(Md,0_d) = 0, V'{AIdAd) = and 
4>Di < then there exists a NPSW solution of the energy integral (fT7| for at least 
one value of M > AId. Hence in the later case double layer solution is unable to 
restrict the occurrence of NPSWs for AI > AId. From this consideration, it is also 
clear that the occurrence of NPSWs is restricted by Ale < AI < AIjj, provided that 
there exists one and only one double root of the equation V{AI]j,(j)) = for the 
unknown 7^ 0, otherwise NPSWs can exist for all AI > AId. 

Therefore, the double layer solution of the energy integral (jl7|) plays an important 
role to determine the upper bound of the Mach number for existence of either PPSWs 
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or NPSWs. For the present problem, PPSW exists if Mc < M < M^ax but from 
the discussion of section 15.11 it is not clear whether the energy integral (|17p provides 
a PPDL solution for some M = such that < M4 < M^ax- Similarly, 
from the discussion of 15.21 it is not clear whether the energy integral ((T7)) provides 
a NPDL solution for some M = M^. So, in the next subsection, we shall analytically 
investigate the existence of double layer solution of the energy integral (fTT]) . 

5.3. Analytical study of the Double layer solution of the Energy Integral 

From the condition (b) as given in section |3] for the existence of double layer solution 
of the energy integral P7|) . we must have a non-zero (0 7^ 0) such that following 
conditions are simultaneously satisfied: 

y(0)=O,r((/)) = O,T/"(0)<O. (25) 

Using - (pOj) . the first equation, the second equation and the third inequality of 
(|25p can be written, respectively, as 

V{(t>) = (1 - nj)Af2 + a- n^iia - 20 - 2anj) 

-5 = 0, 



V'{(f) = n. 



dS 



0, 



where 



Eliminating Ui from ([26| and (|27|) . we get 

, dS 



(26) 
(27) 

(28) 
(29) 



—a. 

It can be easily checked that 4 
(p, (|30p can be written as 

M2 = hicj,), 

where 



dS 



(30) 

if and only if Ue = 1, and consequently, for non-zero 

(31) 



dS 



Using dSI]), from ^ and 



3a - 20- 2a (^f) 



fl{l ~ rie) 
we, respectively, get 

^/2h{^ dS 



= 0, 



where 



(32) 

(33) 
(34) 

(35) 
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Now suppose (f> — (f>4i ^ solves (1551) . then the double layer solution of the energy 
integral PT)) exists at M = M^/ = y'hjjpdi) having amplitude {(pdil if the following 
conditions are satisfied: 

h{cj)di) - Ml > 0, (36) 
g-iM-(f>di>0, (37) 

To derive the condition (|57)) . we have used the following restriction on (j> : (j) < = 
i(M — ^/iaf' . Again the condition states that the function ri{(j)) decreases in a 
neighbourhood of (fidi and r]{(j)di) = 0. Therefore if 0^; > 0, 77(0) changes sign from 
positive to negative when it crosses the point (j)di from left to right. However, for 
4>di < 0, rj{(f>) changes sign from positive to negative when it crosses the point (pdi from 
right to left. In any case, ry(0) vanishes at (f> — (pdi- 



5.4- Algorithm to find negative (positive) potential double layer solution 

From this theoretical discussion, it is simple to make a numerical scheme to test 
whether the energy integral provides a double layer solution. The simple algorithm 
for the existence of double layer solution can be written as follows. 
Step - 1: Find ft.((/)) defined by ([211) for </) < (0 > 0). If h{,p) is negative for all 
< (0 > 0) then the system does not support negative (positive) potential double 
layer. Else follow the next step. 

Step - 2: If y/g+{<j)) — \J 5-(0) — i> be real for < (0 > 0) then go to next step. 
Otherwise, the system does not support negative (positive) potential double layer. 
Step - 3: Set up a numerical scheme to find all possible real negative (positive) roots 
of of ((33)) for all admissible values of the parameters. To find these roots of at 
some fixed values of the parameters, let free and take all those 0's where the function 
77(0) changes its sign from positive to negative when it crosses the point (\) from right 
to left (left to right) and denote the set of all these (/)'s as 2?_ (2?+)- Obviously, there 
may exist other roots of of (I33p for unknown such that 77(0) changes its sign from 
negative to positive, however, the theoretical discussions suggest that these roots are 
beyond the scope of the present numerical scheme. 
Step - 4: Using ([HI]) find the Mach numbers M at each of V- {!)+). 
Step - 5: Find the largest (smallest) value (pn {4>d) of </> from I?_ (2?+) such that 
conditions ([551 - (1571) hold good. Then, obviously (po {4>d) is the amplitude of NPDL 
(PPDL) solution. Next use (|5T|) to obtain the Mach number Mu {Md) corresponding 
to NPDL (PPDL) solution having amplitude |0d| {(pd)- 

With the help of this algorithm, we want to investigate numerically the following 
facts for the present system: 

(i) Whether the present system supports PPDL and/or NPDL solutions, (ii) Whether 
the double layer solution, if exists, can restrict the occurrence of all solitons of same 
polarity, i.e., whether the double layer solution is the ultimate solution of all solitons 
of same polarity of the present system. In this connection, it is important to remember 
that for any double layer solution, there must exists at least one sequence of solitons 
of same polarity converging to the double layer solution, i.e., the amplitude of any 
double layer solution acts as an exact upper bound of the amplitudes of at least one 
sequence of solitary waves of same polarity. 
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Table 1. The values of (po and the corresponding values Mo calculated from 
II31I I have been tabulated for different values of o, ^ and which have been used 
in figure [2l 





a 




/3i 


Md 


<t>D 


Vi: 


0.9 


0.2 


0.6 


5.10476181 


-1.9397924 


V2: 


0.9 


0.5 


0.4 


2.49995351 


-1.5418283 




0.9 


0.5 


0.6 


3.25924103 


-1.5492367 



To investigate the existence of PPDL, r]{<j)) is plotted against <j) for > in figure[T] 
for different values of /i and /3i. Figure [IJa) and figure [rfb) show that ri{(j)) is either 
remain positive throughout (/) > or changes sign from negative to positive when it 
crosses positive (j) axis from left to right. In figure[IJa), corresponding to /3i = 0.4 
changes sign from negative to positive when it crosses positive (p axis from left to right. 
Similar facts have been observed in figure[lja) for ri((j)) corresponding to /3i = 0.6 and 
in figure [TJb) for 7/(0) corresponding to /3i = 0.6. Consequently, there does not exist 
any PPDL solution. For any admissible values of the parameters involved in the 
system, it can be easily verified that the system does not support any PPDL solution 
with the help of plotting ri{(f)) against (j) for > 0. Since the system does not support 
any PPDL solution, we can conclude that Mmax is the only upper bound of M for the 
occurrence of PPSWs. Thus PPSWs exist whenever Mc < M < Mmax- 

To investigate the existence of NPDL, 77(0) is plotted against <j) for (p < Q m 
figure [2] In figure [2fa), 77(0) corresponding to /3i — 0.2 and /3i — 0.4 remain positive 
throughout <j) < and thus there does not exist NPDL solution for f3i — 0.2 and also for 
/3i — 0.4 with a ~ 0.9, /i — 0.2. It can be easily checked that there does not exist any 
NPDL solution for < /3i < 0.4 with a — 0.9, /i = 0.2 by simply drawing 77(0) against 
(j). However ri{(f>) corresponding to /3i = 0.6 changes sign from positive to negative 
when it crosses negative cj) axis from right to left and fulfill all the conditions of our 
algorithm for the existence of NPDL solution. A similar interpretation of figure [Hb) 
shows that NPDLs exist for /3i > 0.4 with a — 0.9, /i = 0.5. Thus for larger values 
of 1-1, NPDL solution is expected at lower values of /3i. Therefore, for denser electrons 
in the background plasma, we require less amount of fast energetic electrons to get 
NPDL. In other words, electron density depletion restrict the occurrence of NPDL. 

Still it is not clear whether the existence of NPDL can restrict all NPSWs of 
the present system, or in other words, whether Md is the upper limit of M for the 
existence of all NPSWs of the present system. From figure [21 one can find three set 
of values of the parameters a, fi and /3i such that NPDL exist. It is easy to find 
at which r/(0) changes sign from positive to negative when it crosses negative axis 
from right to left and using (1311) we can find three values of Mu corresponding to 
three different values oi (j)£) for three set of values of the parameters a, ji and /3i as 
shown in figure [21 For clarity we have tabulate these values in table [H Our algorithm 
will be verified if we can confirm the occurrence of NPDL solutions at those values of 
parameters. 

In figure [3l[a), V{4>) is plotted against (j) for three set of values (denoted as Vi, V2 
and V3 in table [T]) of a, /i, /3i and Md- Our aim is to show the amplitudes obtained 
from figure [3l[a) are exactly the same as obtained from table [H Each of the curves 
in figure [3ja) shows the existence of a NPDL solution. Moreover, the amplitudes of 
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these double layers are exactly the same as obtained in table [T] Hence our algorithm 
regarding the double layer solution is correct. From table [U we have some typical 
observations regarding the amplitude of NPDLs. For fixed a and the amplitude of 
double layer increases with f3i. Again for fixed values of a and /3i, the amplitude of 
double layer decreases with increasing /i. In other words, NPDL gets stronger with 
fast energetic electrons. However, for any fixed non-zero value of /3i and for any value 
of a, one can get stronger NPDL by adding more electrons on the dust grain surface. 

Now we are in a position to investigate whether the NPDL solution can restrict 
the occurrence of all NPSWs of the present system. For this purpose, we explore 
figure [Slja) beyond (j)D and obtain figure [^b). In this figure, we have drawn the 
same set of curves as in figure [3ja) with an exception that we have extended the 
range of 4> axis far away from (p = 0. We see that after making a double root at 
(j) = (j)D: there exists an (poi < (I'd such that V{Mo, i'Di) = 0. Thus according to our 
theoretical discussions in section (HUl there exists a M > Md such that NPSW exists. 
Therefore, the NPDL solution cannot restrict the occurrence all NPSWs of the present 
system and consequently, Md cannot act as an upper bound of M for the existence 
of all NPSWs of the present system. Actually, the present system supports very large 
amplitude NPSW for aU M > Md- To justify this fact, in figure ID^a) and figure IDJb), 
V{(j)) is plotted against (p for three different values of M, viz., Md, Md — 0.005 and 
Md + 0.005. Figure |3i;b) shows that at M = Md, there exists a NPDL of amplitude 
\(j)D\ whereas, for M — Md — 0.005, there exists a NPSW of amplitude less than |0d|. 
However, figure Hfb) shows that the equation V{Md + 0.005, (j)) = has no real root 
of 4> in the neighborhood oi (p = (j>D- From figure lU^a), we see that y(M, again 
vanishes at cj) = (pi, (j> — (j>2 and (p — (p^, respectively, for M = Md — 0.005, M — Md, 
and M = Md + 0.005. However, the roots (p = (pi and = 02 of V{M, 0) = 
corresponding to M = Md — 0.005 and M = Md are unable to give any solitary wave 
solution, whereas the root = 03 of ^(Mc + 0.005, </>) = gives a NPSW of amplitude 
much greater than that of NPDL at M = Md as well as NPSW at M = Md - 0.005. 
Thus there is a finite jump in amplitudes between two NPSWs at M — Md ~ 0.005 
and at M = Md + 0.005 separated by the NPDL at M — Md- This is not a new result, 
the same result has also been observed in some recent works [JU with different 
plasma environments. Mathematically, it is simple to prove the following property: 

Property: If there exists two types of NPSWs (PPSWs) separated by a NPDL 
(PPDL) then there is a finite jump between the amplitudes of two types of NPSWs 
only when |^ < for all M > and for all 0< O (0 > O). 

For the present problem, it is easy to check that 



for all M > and for all (j> ^ 0. Thus all the conditions of the property are satisfied 
but in the positive potential side, there does not exist any jump in amplitudes between 
two solitary waves. More specifically, we have not found any PPDL solution which 
separates two types of PPSWs. 

In figure Hfc), profiles of NPDL has been shown at M — Md- In figure lUJd), 
profiles of NPSWs have been shown at M ^ Md - 0.005 and at M = Md + 0.005, 
respectively. The profile in figure ^c) corresponding to M = AId is an usual 
double layer profile. However the solitary wave profile in figure H^d) corresponding to 
M — Afp + 0.005 is an unusual one; its like a dais-type solitary wave profile. The jump 




(39) 
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between the amplitudes of two NPSWs separated by the NPDL is much prominent 
here. 

In the above discussions, we have demonstrated the possible existence of solitary 
structures for some particular values of the parameters of the problem without making 
any delimitation of the compositional parameter space for the existence of such 
nonlinear structures and consequently, we are unable to produce complete scenario 
of the present problem. So, it is desirable to construct the entire solution space 
or compositional parameter space showing the nature of different solitary structures 
present in the system. In the next section, we have considered different solution spaces 
of the energy integral p7| with respect to Pi. 

6. Different solution spaces of the Energy integral 

Figure H] - figure [S] are the different compositional parameter spaces with respect to 
Pi showing nature of solitary structures and all these figures are aimed to show the 
solution spaces of the energy integral pT| with respect to Pi . To interpret figure |6] 
- figure [9l we have made a general description as follows: solitary structures start to 
exist just above the lower curve M = Mc- For any admissible range of the parameters 
there always exists at least one M > Mc such that NPSW exists thereat. Mmax 
is the upper bound of M for the existence of PPSWs, i.e., there does not exist any 
PPSW if M > Mjnax- More explicitly, if we pick a Pi and goes vertically upwards, 
then all intermediate M bounded by M = Mc and M = Mmax would give PPSWs. 
The curve M = M^ax also restrict the coexistence of both NPSWs and PPSWs, 
however the curve M = M^ax is unable to restrict the occurrence of all NPSWs of 
the present system, i.e., there exists NPSW for all M > Mmax- At any point on the 
curve M — Mo there exists a NPDL solution. But this NPDL solution is unable to 
restrict the occurrence of all NPSWs of the present system. As a result, we get two 
different types of NPSWs separated by the NPDL solution, in which occurrence of 
first type of NPSW is restricted by Mc < M < Md whereas the second type NPSW 
exists for all M > Md- We have also observed a finite jump between the amplitudes 
of NPSWs at M ^ Md - ei and at M = Md + €2, where Q < ti < Md - Mc and 
€2 > 0, i.e., there is a finite jump in amplitudes of the NPSWs above and below the 
curve M = Md- Now we want to define the cut off values of ^ and Pi, which are 
responsible to delimit the solution space. 

Pic '• Pic is a cut-off value of Pi such that NPDL starts to exist whenever Pi > Pic for 
any value of /i lies within the interval < /i < /iy, i.e.. Pi = Pic is the lower bound 
of Pi for the existence of NPDL solution. Thus, Pic is the minimum proportion 
of fast energetic electrons such that maximum potential difference occurs in the 
system and the value of Pic depends on the number of electrons residing on dust 
grain surface. 

(Up : /ip is a cut of value of fi such that Mmax does not exist for any admissible value 
of Pi if fi lies within the interval < 11 < fip, i.e., ii fJ- > fJ-p, there exists a value 
PI of Pi such that Mmax exists at Pi — PI, moreover, if PI > 0, then Mmax exists 
for all Pi lies within the interval < Pi < Pl- 

Pia '• Pia is a cut-off valuc of Pi such that Mmax exists for all < Pi < Pia whenever 
/X > fip. Consequently, Pi = Pia is the upper bound of Pi for the existence of 
PPSW. 
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Now, if Pic > Pia, then there exists an interval /3ia < /3i < /?ic in which neither 
Mjy nor M.fnax exist and consequently, we can define cut-off values fiq and fir of fi as 
follows: 

Hq : iiq is another cut-off value of /i such that for all fJ-p < fJ- < iJ-q, neither M^ax 
nor Md exist whenever /3iq < /?i < /3ic, i.e., for all fJ^p < i-t < fiq and for all 
Pia < I3i < file only NPSWs exist for aU M > M^. 

another cut-off value of fi such that for all /x,. < /i < ^t, the curve 
M = Md tends to intersect the curve M = M^. at the point Pi = /?ic. 

From the definition of fip, fiq and /i^, we can numerically find the values of /ip, fiq and 
/ir for any value of a. The numerical solution is shown graphically in figure [5l From 
this figure we see that for any value of a, we can partition the entire interval of /x in 
the following four subintervals: (i) < /i < /ip, (ii) fip < n < Hq, (iii) fiq < H < Hr and 
(iv) fir < 1^ l£ f-T- In these subintervals of fi, we have qualitatively different solution 
space of the energy integral ([T7)) with respect to /3i. The solution spaces have been 
shown through figure [6]- figure [9] for four different subintervals of fi. 

Before going to discuss the solution spaces in details, the variations of the curves 

Pi = Pia ( ) and Pi = Pic{ ) have been plotted against fi in figure fTOl to 

demonstrate the solution spaces in true physical sense. In this figure, actually we 
have shown those two curves, viz., Pi — Pia and Pi = Pic which are responsible to 
divide fi into several subintervals. A closer look of the figure [5] - figure IHl suggests 
that figure [10] effectively defined all the solution spaces as shown through figure [6] - 
figure [5] in more compact form provided that we have sound knowledge regarding the 
appropriate bounds of the Mach number for the occurrence of different types (nature) 
of solitary structures of the present system. Using the theory as presented in section[51 
one can easily set a numerical scheme to find the appropriate bounds for the occurrence 
of different types (nature) of solitary structures of the present system. In figure [TOl we 
have used the following terminology. C-N: region of coexistence of both PPSWs and 
NPSWs for Mc <M < Mmax and only NPSWs whenever M > M^ax] C-N-D: region 
of coexistence of both PPSWs and NPSWs for Mc < M < M^ax and only NPSWs 
whenever M > Mmax with a NPDL at some M ^ Md > Mc] N-D: region of existence 
of only NPSWs whenever M > Mc with a NPDL at some M ^ Md > Mc. We have 
found fip, fiq and fir lies in the neighborhood of ^ = 0.14, fi = 0.36 and ft = 0.44, 
respectively for a = 0.9. From figure [TUl we have the following observations. 

For Pi = 0, i.e., for isothermal electrons, for < /x < /ip, only NPSWs exist for all 
M > Mc and coexistence of both NPSWs and PPSWs are possible for fip < fi < fiT 
whenever Mc < M < Mmax, whereas NPSWs exist for all M > Mmax- In presence of 
isothermal electrons the system does not support any double layer solution. Similar 
facts can also be observed by considering figure [HI - figure [S] at the point Pi = 0. 
So, if fi exceeds the critical value fip, PPSW starts to exist for Mc < M < Mmax 
and attains its maximum amplitude at M = Mmax but even for increasing fi for 
fip < fi < fiT, the PPSW can not acquire enough strength to make a PPDL even 
at M = Mmax- Actually, From the charge neutrality condition ([7]), we see that n^o 
is a constant. Consequently, we cannot inject positive charge from outside or we 
cannot increase the equilibrium ion number density n^o and this is the reason that 
PPSW cannot acquire enough strength to make a PPDL even at M = Mmax - On the 
other hand, we can increase or decrease the quantities Z^, rido, rico in such way that 
the charge neutrality condition (O holds good for constant equilibrium ion number 
density riio. But in any case, the amplitude of the NPSW steadily increasing for 
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increasing Mach number M > M^. Actually, we are unable to restrict the occurrence 
of NPSW for the present system, i.e., we have not found any upper bound of the Mach 
number which can restrict the occurrence of NPSW and this is the reason that NPSW 
cannot make a NPDL at any point of the compositional parameter space. However, 
to discuss the formation of double layer from physical point of view, we consider the 
following simple mathematics. Suppose (pp is the amplitude of PPSW at any point of 
the compositional parameter space, where we have used the following terminology: if 
there does not exist any PPSW at some point of the compositional parameter space, 
then (j)p — 0. Therefore, 4>p is well defined as the amplitude of PPSW at any point of 
the compositional parameter space. Similarly, one can define ^at, as the amplitude of 
NPSW at any point of the compositional parameter space, i.e., if there does not exist 
any NPSW at some point of the compositional parameter space, then c/jn = 0. From 
simple mathematics, we get 



From inequality (^01) . it is clear that one can get a PPDL solution at a point of the 
compositional parameter space if | (/)p | — | (/)Ar | is maximum with \(t>p\ > \(t>N\ whereas one 
can get a NPDL solution at a point of the compositional parameter space if |(/)Ar| — |(/)p| 
is maximum with \(f>N\ > {(f'pl- For isothermal electrons, it can be easily checked that 
the potential difference (|(/)p| — \4in\ with \(j>p\ > \4in\) for the formation of PPDL can 
not attain any maximum value at any point of the compositional parameter space. 
Similarly, the potential difference (|0Ar| — |0p| with \4>n\ > \4'p\) for the formation of 
NPDL can not attain any maximum value at any point of the compositional parameter 
space. So, for isothermal electrons, the present system does not support any double 
layer solutions. 

For non-zero /3i , the solution space as obtained in figure [10] can be partitioned 
as follows: (i) < < fip-. For < /3i < /3ic, only NPSWs are possible for all 
M > Mc, whereas for fHic < Pi < Pit, NPSWs are possible for all M > Mc except 
M = Md{> Mc). (ii) ^J.p < fJ- < fJ-q-- For < < /3ia, coexistence of both NPSWs 
and PPSWs are possible whenever Mc < M < M^ax and only NPSWs exist for all 
M > Mmax- For Pia < Pi < Pic, Only NPSWs are possible for all M > Mc- For 
Pic < Pi < Pit, NPSWs are possible for all M > Mc except M = Md(> Mc). 
(iii) iJ,q < n < fJ,r- For < /3i < Pic, coexistence of both NPSWs and PPSWs are 
possible whenever Mc < M < Mmax and only NPSWs exist for all M > Mmax- For 
Pic < Pi < Pia, coexistence of both NPSWs and PPSWs are possible whenever 
Mc < M < Mmax and only NPSWs exist for all M > Mmax except the point 
M = Md{> Mmax)- For Pia < Pi < PiT, NPSWs are possible for all M > Mc 
except M = Afp)(> Mc)- (iv) /ir < M < f^T- For < Pi < Pic, coexistence of both 
NPSWs and PPSWs are possible whenever Mc < M < Mmax and only NPSWs exist 
for all M > Mmax- For Pic < Pi < Pu, coexistence of both NPSWs and PPSWs are 
possible whenever Mc < M < Mmax and only NPSWs exist for all M > Mmax except 
M = Mb{Mc < Mb < Mmax)- For Pu < Pi < Pia, coexistence of both NPSWs 
and PPSWs are possible whenever Mc < M < Mmax and only NPSWs exist for all 
M > Mmax except the point M = Md{> Mmax)- For Pia < Pi < Pit, NPSWs 
are possible for all M > Mc except M = Md{> Mc)- In all these solution spaces, 
whenever Md > Mc or Mp, > Mmax, at the point M = Mp,, one can always find a 



(j)p ^ (j)N = \4>P - 4>n\ > |0p| - |0Af| 




(40) 
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NPDL, whereas one can find the coexistence of a NPDL and a PPSW at the point 
M = Md whenever < Md < Mraax- 

Therefore, from the above discussions, it is clear that for any physically admissible 
value of /X, i.e., < /j, < /zt, there exists a non-zero value of /3i such that the present 
system supports a NPDL solution for some M — Md > Mc- So, again from charge 
neutrality condition Q , we see that if the density of electrons increases up to a certain 
value Hr (see figure [TO|l , minimum energetic electrons (small value of /3i ) can produce 
NPDL whereas if the density of electrons tends to zero (almost depletion of electrons) , 
more energetic electrons (higher value of /3i) are required to form a NPDL solution. 
So, we see that /3ic exists for any physically admissible value of /i, i.e., < /i < /it and 
consequently, the present system supports NPDL solution for some Af = AId > Ale if 
Pic < /3i < Pit- 

Again, there does not exist any f3ia for < /i < /ip and consequently, coexistence 
of both NPSWs and PPSWs are not possible even when nonthermal distribution of 
electrons becomes isothermal one, i.e., when /3i = 0. It is also important to remember 
that if the value of /3i increases, negative potential is stronger than positive potential 
and consequently, instead of getting PPSW, from inequality (|40| . one can get a NPDL 
solution. 

For < /I < fJ-T, Pia always exists and increases with increasing ^. Consequently, 
for this interval of /i, solitary structures of both polarities exist provided that 
< ^1 < pia and the region of coexistence of both NPSWs and PPSWs with respect 
to the nonthermal parameter increases with increasing /i lying within /ip < /i < /^t- 
Moreover, for the values of /i lying within /Xp < ^ < /^t, in. fact, fJ-r l£ fJ- l£ /^t, if 
< /3ic < /3ia, the NPSW is stronger than PPSW for < ;3i < Pic and at Pi = Pic, 
inequality pOj) holds good for > \(j)p\ and consequently, we have not only get 
a NPDL solution a. Pi = Pic but also get a weaker PPSW at the same point of the 
compositional parameter space when M = Mjj. 

For Hp < jjL < jjbq, /3ic is always greater than /Sia- Consequently, for /3ia < Pi < 
Pic, only NPSWs exist for all M > Mc whereas for jiq < ^ < ^t, Pic is always less 
than /3ia and all types of solitary structures are possible for the present system. All the 
facts can also be verified by considering figure H]- figure 13 The physical interpretation 
for the formation of solitary structures in this case can be demonstrated through 
charge neutrality condition Q, inequality (j40|) and either considering the figure [TOl or 
more explicitly the figures figure |6] - figure HI 

Therefore, this figure[TO]is actually the graphical presentation of different solitary 
structures with respect to different subintervals of /i within the admissible interval of 
the nonthermal parameter /?i . 

Finally, considering any solution space, we can get new results and physical ideas 
for the formation of solitary structures if we move in the solution space along the family 
of curves parallel to the curve M = Mc- For example, we shall consider the solution 
space with respect to the nonthermal parameter /3i for /i^ < < Mt(< 1) and if we 
move in the solution space along the family of curves parallel to the curve M = Mc, it 
is simple to understand the mathematics as well as physics for the formation of double 
layer solution and it is also simple to understand the relation between solitons and 
double layers. To be more specific, solution space for the present system with respect 
to Pi for fir < /i < IJ-t{< 1) has been presented in Fig. [TlTa'). in which the curve 
M = Mc is omitted from the solution space as presented in figure [HI Now consider the 
family of curves parallel to M = Mc- For instance, consider one such parallel curve 
for M = Mc + 0.04 as shown in Fig. [TlTa). In this figure, Pp is the value of Pi where 
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the curve M = Mjj intersects the curve M = A/^ + 0.04, whereas /Sg is the value of /3i 
where the curve M = M^ax intersects the curve M = Mc + 0.04. Fig. [TIT a) can be 
interpreted in the same way of figure[S]with /3ic replaced by /3p and replaced by /3g. 
However, in Fig. [TIT a) the solitary structures exist along the curve M = Mc + 0.04, 
specifically, (i) both NPSW and PPSW coexist for < /3i < /3<„ (ii) only NPSW 
exists for < /3i < I3it and (iii) at = ^p, a PPSW coexists with a NPDL. The 
variation of amplitude of those solitary waves along the curve M = Mc + 0.04 for 
< Pi < Pit have been shown in Fig. ITlT b). This figure shows that the amplitude 
of NPSW decrease with increasing /3i for /3p < /3i < and the amplitudes of NPSWs 
are bounded by the amplitude of NPDL at (ix — (3p. Again, the amplitude of PPSW 
increases with increasing /3i for /3i > /3p having minimum amplitude at /3i = /3p. 
Moreover, the Fig. [TIT b) shows that along the curve M = + 0.04, the amplitude of 
NPSW increases with decreasing M along the curve M = Mc + 0.04 for [3p < [3i < (3q 
and ultimately, these NPSWs end with a NPDL at = w 0.3975. Therefore, 
the solitons and double layer are not two distinct nonlinear structures, i.e., double 
layer solution, if exists, must be the limiting structure of at least one sequence of 
solitons of same polarity. More specifically, existence of double layer solution implies 
that there must exists at least one sequence of solitary waves of same polarity having 
monotonically increasing amplitude converging to the double layer solution, i.e., the 
amplitude of the double layer solution acts as an exact upper bound or Least Upper 
Bound (lub) of the amplitudes of the sequence of solitary waves of same polarity. 
However, we have seen in the literature that when all the parameters involved in 
the system assume fixed values in their respective physically admissible range, the 
amplitude of solitary wave increases with increasing M and these solitary waves end 
with a double layer of same polarity, if exists. Here it is important to note that M 
is not a function of the parameters involved in the system but is restricted by the 
inequality Mc < M < Mn, where M = Mjj corresponds to a double layer solution. 
So we cannot compare this case with the case of M = Mc + 0.04, since Mc is a 
function of the parameters involved in the system and consequently, monotonicity of 
Mc entirely depends on a parameter when the other parameters assume fixed values 
in their respective physically admissible range. But the solitons and double layer are 
not two distinct nonlinear structures. Therefore, double layer solution, if exists, must 
be the limiting structure of at least one sequence of solitons of same polarity. In Fig. 
[TlT b). by the vertical line with both sided arrow, we mean, the amplitude of NPDL 
at /3i = /3p for M = Mc + 0.04. For a = 0.9 and fi — 0.5, the critical values /3p and /?, 
lies in the neighborhood of f3i = 0.3975 and Pi = 0.4623, respectively. Now, for the 
formation of NPDL solution on the curve M — Mc + 0.04 at /3i — /3p, the negative 
potential (absolute value) must dominate the positive potential in the neighborhood of 
the point /3i = Pp and the potential difference (with respect to negative potential) must 
be maximum thereat. From Fig. [TlT b). it is clear that the negative potential (absolute 
value) dominates the positive potential in a right neighborhood of /?i = /?p and the 
potential difference (|0Ar| — |0p| along with \(I)n\ > \(t>p\) is maximum thereat. This 
figure shows the existence of NPDL solution at Pi = Pp, which is already confirmed 
in Fig. [TTT a). More specifically, from the inequality (j40|). it is clear that one can 
get a PPDL solution at a point of the compositional parameter space if \4>p\ — \(I>n\ 
is maximum with \(j)p\ > |07v| whereas one can get a NPDL solution at a point of 
the compositional parameter space if \4>n\ — |0p| is maximum with |(/)jv| > \(j)p\. 
From figure ITTT b) . we have found that \4>n\ — \4'p\ is maximum with \(j)N\ > \4'p\ 
Pi = Pp, and consequently, we can have a NPDL solution at /3i = /3p w 0.3975 for 
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M = Mc + 0.04. Next we consider the curve M — M^. + 0.08 parallel to the curve 
M = Mf, as shown in Fig. [T^ a). Here also, /3p and /3q are defined in the same way as 
in Fig. [TlTa). However, from Fig. [T^ a). we see that there does not exist any PPSW 
along the curve M = Mc + 0.08 for /3p < /3i < /3it- So according to the terminology 
(j)p = along the curve M = Mc + 0.08 for /3p < /3i < /3it- Consequently, from 
inequality (gO]), we have NPDL solution at /3i = /3p w 0.4234 for M = Mc + 0.08. This 
fact is clear from figure [T2l[b) . 

7. Summary and Discussions 

In the present paper, we have investigated DIASWs and DIADLs in a dusty plasma 
system consisting of adiabatic ions, nonthermal electrons and negatively charged dust 
grains. Investigations have been made by going through the entire solution space of 
the energy integral by considering the entire range of the parameters involved in the 
system. Our aim is to delimit the parameter /3i depending on the nature of existence 
of DIASWs and DIADLs. 

Therefore, for any physically admissible values of the parameters of the system, 
specifically, for any value of fi and any value of NPSW exists for all M > Mc 
except M = AId, where Mc is the lower bound of Mach number M, i.e., solitary 
wave and/or double layer solutions of the energy integral start to exist for M > Mc 
and Md is the Mach number corresponding to a NPDL solution. However, if the 
parameter fi exceeds a critical value /Xp, PPSWs exist for all /ip < /i < 1 whenever the 
Mach number lies within the interval Mc < M < M^ax^ where M^ax is the upper 
bound of M which is well-defined only when the system supports PPSWs. Therefore, 
the coexistence of both PPSWs and NPSWs is possible for all /ip < /i < 1 whenever 
Mc < M < M,nax: but NPSWs still exist for M > M^ax- 

For nonthermal electrons, NPDL starts to occur whenever the nonthermal 
parameter exceeds a critical value. However this double layer solution is unable to 
restrict the occurrence of NPSWs. As a result, two different types of NPSWs have been 
observed, in which occurrence of first type of NPSW is restricted by Mc < M < Mjj 
whereas the second type NPSW exists for all M > Mu, where Md is the Mach 
number corresponding to a NPDL. A finite jump between the amplitudes of NPSWs 
at M = M^) — e and at M = Mjj + e has been observed, where e is a sufficiently 
small positive quantity. The amplitude of NPSW for M > Mjj is much greater than 
the amplitude of the NPDL solution at M = Mjj as well as the amplitude of NPSW 
for M < Md, i.e., there is a jump in amplitudes of the NPSWs above and below the 
curve M = Mo- However, there is no jump in amplitudes of NPSWs above and below 
the curve M = M^ax- 

In most of the earlier works, dust ion acoustic solitary structures have been 
investigated with the help of Maxwellian velocity distribution function for electrons. 
However, the dusty plasma with nonthermally/suprathermally distributed electrons 
observed in a number of heliospheric environments [131 HH HOI IHl HH US] • Therefore, 
the present paper gives the complete scenario of dust ion acoustic solitary structures in 
a dusty plasma system in which lighter species (electrons) is nonthermally distributed. 
This paper is also helpful for understanding the formation of dust ion acoustic solitary 
structures from different physical and mathematical aspects. 
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Figure 1. r]{(p) is plotted against cp for a = 0.9 with /i = 0.2 in (a) and fi = 0.5 
in (b). In each of (a) and (b), three curves have been drawn corresponding to 
three different values of /3i, viz., /9i are 0.2[ ], 0.4[ — ] and 0.6[— ■ — ]. 




Figure 2. r]{(p) is plotted against cp for a = 0.9 with = 0.2 in (a) and ix = 0.5 
in (b). In each of (a) and (b), three curves have been drawn corresponding to 
three different values of viz., are 0.2[ ], 0.4[ — ] and 0.6[— ■ — ]. 
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Figure 3. V{4>) is plotted against </> for three set of values Vi , V2 and V3 presented 
in table[Tl In figure (a) all thie curves (Vi , V2 and V3) show the existence of NPSW 
solutions. In figure (b) the same curves have been plotted in an interval of </> close 
to </> = 0. 
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Figure 4. V{cp) is plotted against cp for three different values of M, viz., Mo[ — 

], Mo - 0.005[ ] and Mp + 0.005[- • -] in (a) and (b). In (b) we have 

shown the region of cp from </> = —2 to = 0, whereas in (a) a region of 4> from 
cp = —13.5 to (/> = —12.5 has been shown. The NPDL profile corresponding to 
M = M]j has been shown in (c), whereas the profiles of NPSWs corresponding 
to M = Mb - 0.005 and M = Mo + 0.005 have been drawn in (d). Profiles in 
(d) show a finite jump between the amplitudes of solitary wave by going from 
M = Mo - 0.005 to M = Mo + 0.005. In all these four figures we have used 
a = 0.9, iJ, = 0.5 and /9i = 0.4. 
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Figure 5. fip, [iq and arc plotted against a. Prom this figure we see that for 
any value of a we have ^ip < iiq < fXr ■ 
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Figure 6. A graphical presentation of different solitary structures have been given 
with respect to two different subintcrvals of /3i within the admissible interval of 
the Mach number M for a particular value of \x which lies in < /i < /ip. The 

curves M = { ) and M = Af£)( ) are responsible for the occurrence 

of two subintervals of /9i . At any point M = Md one can always find a NPDL 
solution. This solution space has been drawn for a = 0.9 and /j, = 0.1. For 
a = 0.9, we have found Hp lies in the neighborhood of n = 0.14 and /3ic ^ 0.56. 
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Figure 7. A graphical presentation of different solitary structures have been given 
with respect to three different subintcrvals of /3i within the admissible interval of 
the Mach number M for a particular value of fi which lies in fip < jj, < Hq. The 

curves M = Mc ( ), M = Mmax (- • -) and M = Mu{ ) are responsible 

for the occurrence of three subintervals of At any point M = one can 
always find a NPDL solution. The region in < /9i < has been shown in the 
inset. This solution space has been drawn for a = 0.9 and fi = 0.25. For a = 0.9, 
we have found Hp Ri 0.14 and /Xq ~ 0.36 with /3ia ~ 0.253 and /3ic ^ 0.49. 
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Figure 8. A graphical presentation of different solitary structures have been given 
with respect to four different subintervals of /3i within the admissible interval of 
the Mach number M for a particular value of fj, which lies in fiq < ^ < fir- The 

curves M = Mc ( ), M = Mmax (- ■ -) and M = Mu( ) are responsible 

for the occurrence of three subintervals of ^i. At any point M = Mo one can 
always find a NPDL solution. This solution space has been drawn for a = 0.9 and 
H = 0.4. For a = 0.9, we have found Hq » 0.36 and fXr » 0.44 with /9ia « 0.437 
and /3ic « 0.352. 
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Figure 9. A graphical presentation of different solitary structures have been given 
with respect to four different subintervals of Bi within the admissible interval of 
the Mach number M for a particular value of /i which lies in fir ^ fJ- ^ IJ-t- The 

curves M = Mc ( ), M = Mmax (- ■ -) and M = Alo{ ) are responsible 

for the occurrence of four subintervals of ^Si. At any point M = Mo one can 
always find a NPDL solution. This solution space has been drawn for a = 0.9 
and // = 0.5. For a = 0.9, we have found iXr lies in the neighborhood of ix = 0.44, 
and 01a 0.511, ftt fs 0.417, /3ic fs 0.37. 



Dust ion acoustic solitary structures 



27 



a = 0.9 



0.7 



CO. 

■o 0.35 

c 

a 

n 



N-D 


N-D 




N-D ^ ^ 


s 


s 

\ 

\ 

\ 




y"^ C-N-D ^ " ^ 

/ ^ 
/ ^ 

y 

/ 

/ 




N / 


\ 

\ 


/ 

/ 


N 


/c-N 




C-N 



Figure 10. Here fi\ = /3ia ( ) and fi\ = I3\c{ ) have been plotted against 

jx for a = 0.9. In this figure, actually we have shown those two curves, viz., 

I3i = pia ( ) and /3i = /3ic( ) which are responsible to divide into 

several subintervals. This figure is actually the graphical presentation of different 
solitary structures with respect to different subintervals of // within the admissible 
interval of the nonthermal parameter 



Dust ion acoustic solitary structures 



28 




Figure 11. (a) Solution space for the present system with respect to /3i has 
been presented, where the curve M = Mc is omitted from the solution space as 
presented in figure (9] (b) Variation in amplitude (absolute value) of both NPSW 
and PPSW have been shown along the curve M = Mc + 0.04 for fip < fii < /3q. 
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Figure 12. (a) Solution space for the present system with respect to /3i has 
been presented, where the curve M = Mc is omitted from the solution space as 
presented in figure |9] (b)Variation in amplitude (absolute value) of NPSW have 
been shown along the curve M = Mc + 0.08 for /3p < /3i < /3it- 



